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Mirror matter is a promising self-collisional dark matter candidate. Here we study tlie evolution 
of tiiermodynamical quantities in the early Universe for temperatures below ~ 100 MeV in presence 
of a hidden mirror sector with unbroken parity symmetry and with gravitational interactions only. 
This range of temperatures is interesting for primordial nucleosynthesis analyses, therefore we focus 
on the temporal evolution of number of degrees of freedom in both sectors. Numerically solving 
the equations, we obtain the interesting prediction that the effective number of extra-neutrino 
families raises for decreasing temperatures before and after Big Bang nucleosynthesis; this could 
help solving the discrepancy in this number computed at nucleosynthesis and cosmic microwave 
background formation epochs. 
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I. INTRODUCTION 

In 50's Lee and Yang [l] suggested that the existence of 
a hidden sector of particles and interactions with exactly 
the same properties as our visible world (except for right- 
handed, instead of left-handed, weak interactions) could 
restore the total symmetry of physical laws. Many years 
later Foot at al. [4I proposed a model with exact parity 
(mirror) symmetry between the two sectors, that are de- 
scribed by the same Lagrangian and coupling constants, 
and consequently have the same microphysics. Since the 
only interaction linking the two worlds is gravity[23| , mir- 
ror baryons constitute a dark matter candidate in a nat- 
ural way, as they interact with mirror photons, but not 
with ordinary ones. 

The phenomenology of mirror matter was studied in 
several papers (for an extended list see the bibliogra- 
phy of Ref. 3); in particular the implications for Big 
Bang nucleosynthesis 0, Q , primordial structure forma- 
tion and cosmic microwave background [a, lg, l7|,|8|, la 1C|, 
large scale structure of the Universe 
microlensing events (MACHOs) 
tation of DAMA experiment |16l . Il7l | are related to this 
study. 

If the mirror (M) sector exists, then the Universe along 
with the ordinary (O) particles should contain their mir- 
ror partners, but their densities are not the same in both 
sectors. Indeed, the BBN bound on the effective number 
of extra neutrino families implies that the M sector has a 
temperature lower than the O one, as naturally obtained 
in some inflationary models fisj. Then, two sectors have 




different initial conditions, they do not come into ther- 
mal equilibrium at later epoch and they evolve indepen- 
dently, separately conserving their entropies and main- 
taining approximately constant the ratio among their 
temperatures. 

All the differences with respect to the O world can be 
described in terms of only two free parameters, defined 
as: 

1/3 

; p = n'jny,, (1) 

where T (T'), ilf, (ilj,), and s {s') are respectively the O 
(M) photon temperature, cosmological baryon density, 
and entropy density. :24| The bounds on the mirror pa- 
rameters are x < 0.7 and /3 > 1, the first one coming 
from the previously mentioned BBN limit and the sec- 
ond one from the hypothesis that a relevant fraction of 
dark matter is made of M baryons. In general, during 
most of the history of the Universe, we can approximate 
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where now q{T) and q' {T) are respectively the O and M 
entropic degrees of freedom, and the M photon tempera- 
ture T'(T) is a function of the O one. This approximation 
is only valid if the temperatures of two sectors are not 
too different, or alternatively if we are away enough from 
crucial epochs, like e~'"-e~ annihilation But we are 
just investigating the range of temperatures interested 
by this phenomenon, and in fact one of the aims of this 
paper is to exactly compute the trends of these thcrmo- 
dynamical quantities when the previous approximation 
is no longer valid. 

The present study is required for different reasons. 
First of all, we need a detailed study of thermodynam- 
ical evolution of the early Universe in order to obtain 
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a reliable description of the thermal cosmic history in 
presence of M dark matter. Secondly, an accurate study 
of the trend of number of degrees of freedom is neces- 
sary for any future simulation of nucleosynthesis of pri- 
mordial elements in both sectors, that may be compared 
with observations. In particular, it becomes even more 
important in view of the recent proposed interpretation 
of the DAMA/LIBRA experiment in terms of interac- 
tions with M heavy elements [l^, [l7[ ■ Thirdly, claims for 
variations of the effective number of extra neutrino fami- 
lies computed at BBN 1 MeV) and cosmic microwave 
background (CMB) formation (< 1 eV) epochs require 
investigations for possible variations of these numbers as 
consequences of physics beyond the Standard Model (see 
the recent Ref. [l9| and references therein for previous 
works) . 

The plan of the paper is as follows. In next section we 
introduce the thermodynamical equilibrium in the early 
mirror Universe and the equations that govern its evo- 
lution. In section Hill we present the numerical solutions 
of these equations and discuss the results. Finally, our 
main conclusions are summarized in section ITVl 



II. THERMODYNAMICAL EQUILIBRIUM 

We extend the standard theory of thermodynamics in 
the early Universe in order to take into account the ex- 
istence of M particles. 

We consider the Universe as a thermodynamical sys- 
tem composed of different species (O and M electrons, 
photons, neutrinos, nucleons, etc.) which, in the early 
phases, were to a good approximation in thermodynam- 
ical equilibrium, established through rapid interactions, 
in two sectors separately. 

We assume, as usual, that the Universe is homoge- 
neous, and the chemical potentials of all particle species 
A are negligible, i.e. /Lt^ < T [13, HH. The latter as- 
sumption implies the conservation of the total entropy 
of the Universe S. Therefore we can use the equilibrium 
Bose-Einstein or Fermi-Dirac distribution functions and 
calculate the energy density p, the pressure p, and the 
entropy density s for every particle species in thermal 
equilibrium: 
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where gA is the spin-degeneracy factor of the species A, 
the signs -I- and - correspond respectively to fermions and 



bosons, and E = -^/p^ -I- m^, with p the momentum. We 
can use the usual parametrization of the entropy density: 



where 



s{T) = — q{T) 
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fermions 
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is the number of effective entropic degrees of freedom 
(DoF), that group together the DoF for all bosons (gb) 
and fermions (5/). An analogous formalism can be used 
for the total energy density p, which can be parametrized 
as: 



where 



p(T) = -^9iT)T\ 
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is the number of effective energetic degrees of freedom. 



A. Equations 

We calculate the equations which link the O and M sec- 
tor temperatures and thermodynamical quantities; then 
we numerically solve these equations. Once the temper- 
atures are known, it is possible to work out the total 
exact number of DoF in both sectors, which can be, as it 
is common in the literature, expressed in terms of extra- 
neutrino number. We report in mill the results of these 
calculations. 

The presence of the other sector leads in both sectors 
to the same effects of having more particles. As already 
stated, we do not take in account interactions between 
the two sectors besides gravity. This implies that the 
entropies of two sectors are separately conserved, and 
the parameter x is constant during the cosmic evolution: 



1/3 



/ 3\ 1/3 



1/3 



const. ,(10) 



where a is the scale factor. 

As already stressed, O and M sectors have the same mi- 
crophysics; therefore we can assume at first approxima- 
tion that the neutrino decoupling temperature [25j To^ is 
the same in both of them, that is = T^^,. But the 
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temperatures in the O sector when the O and M neutrino 
decouplings take place are different, T^^, ^ Tj^^i, with 
Tdu < Tdu' since a; < 1. Therefore, at T > Tdu {T'j^^), 
O (M) neutrinos are in thermal equilibrium with the O 
(M) plasma and their temperature isTi, =T [Tl — T'), 
while after decoupling it scales as a~^. 

In each sector, shortly after the v decoupling, the e+- 
e~ annihilate because the temperature becomes smaller 
than 2me, which is the threshold for the reaction 7 ^ 
e+e" in both O and M worlds. Thus electrons and 
positrons transfer their entropy to the corresponding pho- 
tons, which become hotter than neutrinos. 

This fact will be used together with the entropy con- 
servations (total and in each sector separately) to find 
the equations that govern the evolution of the M photon 
temperature T' as a function of the O one T. In fact, O 
and M neutrino decouplings are key events, together with 
both e~^-e~' annihilation processes, for the thermodynam- 
ics in this range of temperatures. Once we call Tu^i the 
O world temperature when the M neutrino decoupling 
takes place, we can split the early Universe evolution for 
temperatures T < 100 McV into three phases. 



1. Phase T > T^^' 

Photons, electrons, positrons and neutrinos are in ther- 
mal equilibrium in each sector separately, that is T^, = 
Te=T , = = T'.fH Using equations from © to 
^ for particles in both sectors we are able to calculate 
the DoF number in O or M worlds alone; but to work out 
the total DoF number, that is summed on both worlds, 
we need the M temperature T' as a function of the O one 
T or viceversa. 

As we neglect the entropy exchanges between the sec- 
tors (valid since there are only gravitational interactions 
between them), we can obtain both these functions using 
Def. ([T]) and imposing x = constant in: 



s' + s' + s' 
X = ■ 



97 



, (11) 



where qv — Q and — 2, while qe{T) stands for 



2. Phase Tdu <T < To,y 

At T ~ To^' M neutrinos decouple and soon after M 
electrons and positrons annihilate, transferring their en- 
tropy only to the M photons, and hence raising their 
temperature. In the M sector the entropies of the system 
{e'^ ,y) and of j^' are separately conserved. Nevertheless, 
O photons and neutrinos still have the same temperature 
T. Therefore we have two equations. The first one comes 
from the conservation of entropies in the M sector, so that 
their ratio is equal to the asymptotic value computed at 



high temperatures, when the e 
had not begun yet: 



annihilation process 
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(13) 



and its solution gives Tl as a function of T'. The second 
one is the conservation of ratio of entropies in two sectors: 



[lq,{T') + q^] r^ + lq^Tj- 



(14) 



and its solution, obtained using Eq. (fT3j) . gives T' as a 
function of T. Eq. (fT4|) is the same as Eq. (fTTjl but with 

rp! rpl rpt / rpf 



3. Phase T < To,. 

At T ~ Tdi^ O neutrinos decouple and soon after O elec- 
trons and positrons annihilate. Now also in the O sector 
the entropies of the system (6^,7) and i' are separately 
conserved; therefore we need one more equation to work 
out the O neutrino temperature T^, as a function of the 
O photon one T: 
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with Se{T) defined in (O, and we have used T'^ = = 
T!^ = T', Te = T„ = ^ T. We have also used that, 
since particles and physics are the same in both sectors, 
the equilibrium distribution functions, and thus the spin- 
degeneracy factors, are the same: q^ — q'^, q.y = qiy, q^ = 
ql- 

Equation (llip can be numerically solved in order to 
obtain the function T'(T) for every T > Tdu' ■ 



[lq,{T') + q^] r^ + lq,Tl^ 



'ke{T) 



q^\ T^ + lq.TS 



(17) 



Eq. ([TT]) is the same as Eq. ^ but with ^ = T ^ 

Once both O and M photon temperatures are known, 
it is straightforward to calculate the total energy and en- 
tropy densities; then, reversing equations ([6]) and (O, we 
can work out the entropic {q) and energetic (5) number 
of DoF. Calculations have been made for several different 
values of x, as reported in the following section. 
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FIG. 1; The ratio ^ for several values of x. The asymptotic 
values of this ratio are 1 for both high and low temperatures, 
as expected. 



III. NUMERICAL CALCULATIONS 

The equations we introduced in the previous section 
have been numerically solved. Since x is a free parameter 
in our theory, several values have been used for it - from 
0.1 to 0.7 with step 0.1 or less. In the extreme asymp- 
totic cases T > Td^' ~^otT<$: T<j„„ e± ~ 1 MeV we 
expect q{T) ~ q'[T'{T)]; therefore, using Eq. (O, in this 
limits the ratio of M and O photon temperatures should 
be X, that is ^ ~ 1. Instead, when T^^/ >T> T^nn e± 
we expect q'[T'{T)] < q{T) because the e+-e~ annihi- 
lation takes place before in the M world, leading to a 
decrease of q'(T') and a corresponding increase of T' in 
order to keep constant the entropy densities ratio; thus 
we expect ^ > 1. Later on, when the O electrons and 
positrons annihilate, they make even T increase and thus 
the ratio ^ decreases to the asymptotic value 1. These 

remarks have been numerically verified; the ratio ^ is 
plotted in Fig. [T] for different values of x. 

M e~^-e~ annihilation happens at higher O tempera- 
tures for lower x values, raising before the temperature 
of M photons; this results in a shift of the peaks of Fig.[T] 
towards higher T. In addition, for lower x the difference 
in T between the two annihilation processes is higher, so 
M photons have more time to raise their temperatures 
before O ones start to do the same, and this leads to the 
change of the shape of the curves. 



A. Number of degrees of freedom 

Using equations from ([3]) to ([5]) to reverse (O and ([5]), 
we can work out the total number of O entropic (g) and 
energetic {g) DoF at any temperatures T. We can apply 
the same procedure to work out the standard., that is in 
absence of the M sector, total qstd and gstd^ as well as 



the M ones q' and g', and the influence of a world on the 
other one. 

At first approximation we expect q (g) to have a cu- 
bic (quartic) dependence on x when the temperature 
is not close to the v decoupling and the e^-e~ anni- 
hilation phases: q = qstd{^ + x^), q' = Qstdi^ + x^^), 
g = ffstdCl + a;^)j g' = gstdi^ + x^'^) (for an explanation 
see Ref. Q). 

In the following we present the results of accurate nu- 
merical calculations. In Table U some values are reported 
for special temperatures and several values of x. As ex- 
pected, the total DoF numbers are always higher than 
the standard and increase with x. Moreover, the M sec- 
tor values are higher than the O ones by a factor of order 
x~^ (for q) or x~'^ (for g). 

In Fig. [21^ qstd, 9, are plotted for the intermediate 
value X = 0.5; panel B of the same figure shows the 
corresponding values of g. In the figure q' {g') has been 
scaled by a factor x^ (x*); in this way the asymptotic 
values are the same than the O sector ones because at 
the extremes ^ = x (see Fig. [T]). We can see that g, g, 
q' and g' have similar trends, but, due to T' < T, DoF 
in the M sector begin to decrease before (at higher T). 

In Fig. [3] 5 in O and M sectors are plotted in com- 
parison with the standard for several x values (from 0.1 
to 0.7 with step 0.1). The predicted quartic dependence 
of g on a; at the extremes is proved correct. In panel 
A we note that ordinary DoF with x < 0.3 are practi- 
cally identical to the standard case. In addition, the plot 
shape does not change with x in the O sector, while it 
does in the M one. This sector evolves with temperature 
T' ^ xT < T; therefore, for lower x the number of DoF 
starts before decreasing below the asymptotic value at 
high T. The change of the shape of the plots in panel B 
is related to the same physical processes responsible of 
the analogous effect present in Fig. [l] that is due to the 
e+-e^ annihilation in the M sector. 



B. Number of neutrino families 

We know that the SM contains three neutrino species; 
the possible existence of a fourth neutrino has been in- 
vestigated for a long time, also using BBN constraints. 
This is why in the literature one can often find bounds on 
the number of DoF in terms of effective extra-neutrino 
number AN^, — — 3. In general, the effective num- 
ber of neutrinos N^, is found assuming that all particles 
contributing to the Universe energy density, to the exclu- 
sion of electrons, positrons and photons, are neutrinos; 
in formula that means: 

j(T) = eAT) + + ''-.2N, ■ i^') 

^A-^^ ^'^'-ff -^ ■(^)'. (18) 
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TABLE L Standard and non-standard total DoF numbers for 
several x values in both O and M sectors. Temperatures are 
in MeV. 
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Ni, has been worked out using Eq. p8)) together with 
the results of previous numerical simulations; some data 
are reported in Table [TTl while plots for several x val- 
ues and any temperatures from to 3 MeV are shown 
in Fig. m We stress that the standard values N^, = 3 is 
the same at any temperatures, while a distinctive feature 
of mirror scenario is that the number of neutrinos raises 
for decreasing temperatures. Anyway, this effect is not 
a problem; on the contrary it may be useful since recent 
cosmological data fits give indications for a number of 
neutrinos at recent times higher than at BBN epoch. In 
Ref. [l^ the authors found A^^ = 5.2^2^2 using recent cos- 
mic microwave background (CMB) and large scale struc- 
ture (LSS) data. We can use the aforementioned mirror 
feature together with the results of a previous work on 
(CMB) and (LSS) power spectra [1, [l3|, where author 
studied the dependence of the spectra on the parameters 
X and for a flat Universe with mirror dark matter. 
As can be easily evinced from Figures 11, 12, 14, 15 of 
Ref. , an increase of the effective number of neutrinos 
is well mimicked by an increase of the parameter x, and 
the amounts of the respective increases are in accordance 
with what required to justify the data of CMB and LSS. 
Thus, considering the sum of these effects, i.e. the raise 
in Ni, before and after BBN, and the similarity of CMB 
and LSS spectra of mirror models and standard A'^^ with 
the ones obtained without mirror sector but with larger 
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FIG. 2: Total entropic (panel A) and energetic (panel B) de- 
grees of freedom computed in ordinary and mirror sectors and 
for the standard. The mirror values have been multiplied by 
'^'^ (q) and x"^ (g) to make them comparable with the ordinary 
1. 



ones, since 



T' 
xT 



TABLE II: Effective number of neutrinos in the ordinary sec- 
tor for some special cases. Temperatures are in MeV. 
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Ni,, the mentioned discrepancy naturally disappears. 

Similarly the effective number of neutrinos in the M 
sector can be worked out as: 



N' 



g'eiT') 



^ ■ 2 

8 ^ 



(19) 



Once again these values are higher than the O ones, but 
now by a factor x"^; they have been numerically com- 
puted and some special values are reported in Table IIIII 
For lower x this number can become very high, inducing 
relevant consequences on the primordial nucleosynthesis 
in the M sector, that is highly dependent on this param- 
eter. 
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FIG. 3: Total energetic degrees of freedom in the mirror sec- 
tor computed for several values of x and compared with the 
standard. 



TABLE III: Effective number of neutrinos in the mirror sector 



for some 


special cases. 


Temperatures are in MeV. 
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a. Predictions for special temperatures It is possible 
to work out the asymptotic values of a,tT :s> Tdu and 
T <C T'onn ™ ^ simple way starting from the standard 
values: 

gstd{T > Td.) = 10.75 ; 

gstd{T < Tanned) — 3-36 . 

Without the M sector, we have, as expected: 

10.75 -2-1 -A 



(20) 



iv,(T»r,3,) = 



^ • 2 

8 ^ 

3.36-2 



•2 



While, when the M sector is present, wc can use the pre- 
viously mentioned quadratic approximation g = gstdi^ + 
x'^) at the special temperatures we are considering; hence: 



10.75(1 + a;4) 



= 3 + 6.14 x'^ ; 

3.36(1 + a;-^) - 2 

8 ^ 
~ 3 + 7.40 a;* . 



(22) 



Prom the equations written above we can see that the 
rise AN^ is: 



. 1 
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8 ^ 



±)-7V,(T»Tz3,) 
- 10.75 



3.36 ( ^ 



~ 1.255 • . 



(23) 



= 3. (21) 



This leads to iV,^ higher than the standard in presence of 
the M world and to a further rise of this parameter at low 
temperatures. If we assume a conservative limit on the 
effective number of neutrino families Ni,{T ^ Tb;/) < 4, 
this implies x < 0.635 and N,{T < T<j„„e±) < 4.26. For 
some interesting values of the parameter x we obtain: 

a; = 0.7 ^ N,{T » Tdv) ^ 4.5 and 
7V,(T«T„„„e±)-4.8; 

a; = 0.6 ^ N,{T » Tdu) 3.8 and 
7V,(T<T„„„,±)~3.96; 

a; = 0.3 => Ar^(r > Td^) ~ 3.05 and 

7V,(T«T„„„e±)-3.06. (24) 

These rough estimates are in good agreement with the 
asymptotic numerical values computed before and after 
BBN; nevertheless, the previous detailed analysis of the 
evolution of this quantity is crucial for future studies of 
the nucleosynthesis in both sectors. 



IV. CONCLUSIONS 

In this paper we have investigated the consequences of 
the existence of a mirror sector of particles and interac- 
tions with exact parity symmetry on the thermodynamics 
of the early Universe for temperatures below ^ 100 MeV. 
In particular we studied in detail the evolution of the en- 
tropic and energetic degrees of freedom in both sectors. 

To perform these calculations we first found the sys- 
tem of equations governing the thermodynamical equi- 
libria in the two sectors, and then we numerically solved 
them in order to obtain the mirror temperature T' cor- 
responding to any ordinary temperature T; once these 
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FIG. 4; Effective number of neutrino families in the ordi- 
nary sector for several x values. 



temperatures were known, we were able to work out the 
entropic and energetic densities, and thus the numbers of 
degrees of freedom and the effective numbers of neutrino 
famihes in both sectors. We found that the evolution of 
the ordinary sector is influenced by the mirror particles, 
that provide extra degrees of freedom, and, due to the 
electron-positron annihilation process happening in the 
mirror sector, changes with time for temperatures around 



the BBN epoch. This makes crucial the present investi- 
gation in order to study in details the BBN process in 
presence of mirror dark matter. Furthermore, the special 
feature that models containing the mirror sector change 
the number of equivalent neutrinos during the Universe 
evolution, if considered together with preexisting mod- 
els of cosmic microwave background (CMB) and large 
scale structure power spectra with mirror dark matter, 
furnishes a possible interpretation of the discrepancy in 
the effective number of neutrino families indirectly "ob- 
served" at BBN and CMB epochs. 

In the mirror sector we obtained very high values of 
degrees of freedom, as predicted, and again they are vari- 
able for temperatures below some MeV. These data pro- 
vide the necessary input for a detailed study of BBN 
in the mirror sector. It is required in order to obtain 
the exact amounts of primordial heavy mirror elements, 
that are necessary for a correct interpretation of non- 
linear processes of structure formation and of the new 
DAMA/LIBRA results. 
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